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Abstract 

Neutrino oscillations in the presence of strong gravitational fields are studied. We look at very high 
energy neutrinos (~1 TeV) emanating from Active Galactic Nuclei (AGN). It is observed that spin 
flavor resonant transitions of such neutrinos may occur in the vicinity of AGN due to gravitational 
effects and due to the presence of a large magnetic field (~1 Tesla). A point to note is that matter 
effects (normal MSW transitions) become negligible in comparison to gravitational effects in our 
scenario. 
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1. Introduction 



It is well known that relativistic effects produce a precession of gyroscopes in the vicinity of a source of 
gravitational field It is therefore reasonable to expect that particles will undergo helicity flips in such 
regions. This is particularly interesting in the case of neutrinos because left handed neutrinos traversing a 
strong gravitationaLfield could be converted into unobservable right handed ones, resulting in a decrease 
of the neutrino fluxtl. We will analyze in this paper, following ||, the possibility of encountering resonant 
transitions of high energy Dirac neutrinos produced in Active Galactic Nuclei (AGN), (reminiscent of the 
ones proposed as solutions of the solar neutrino problem 0, |5], ||). The presence of such resonances alters 
the expectation of the neutrino flux on Earth. In our calculations we will consider magnetic, matter and 
gravitational effects. The latter are introduced using the standard couplings of neutrinos to gravity; we 
will not consider possible effects of the modification of the equivalence principle [0 in this paper. 

The search for such high energy neutrinos by the neutrino telescopes (e.g. DUMAND, AMANDA, 
NESTOR, BAIKAL, etc., |§) which are under construction necessitates a clear picture of the expected 
neutrino fluxes from these objects. AGN are by far the strongest sources of ultrahigh energy neutrinos in 
the Universe, producing fluxes expected to be detectable on the earth with present technology ||. These 
objects are believed to be fueled by the gravitational energy of matter accreting onto a supermassive black 
hole (10 4 to 1O 1O M ) at the AGN core, where gravitational energy is converted into luminous energy 
through the acceleration of high energy protons ||. This environment exhibits copious production of 
hadronic and subsequent leptonic high energy byproducts, such as neutrinos. 

The paper is organized as follows. We start with a brief discussion of an AGN model and mechanism 
of high energy neutrino production in it (section 2). In section 3 we study the possibility of neutrino 
transitions between different flavors and spins in the environment of a strong gravitational field using the 
semiclassical approximation. In section 4 we apply the formalism to the interesting case of a Kerr black 
hole; our main results are presented in this section. Finally in section 5 we give our conclusions. Several 
mathematical details are relegated to the appendices. 



2. AGN Model 



Active galactic nuclei have long been recognized as possible sources of high energy signals | 10fl , being 
the most luminous objects in the Universe. They have luminosities ranging from 10 42 to 10 48 ergs/sec, 
corresponding to black hole masses of the order of 10 4 to 1O 1O M , on the natural assumption that they 
are powered by Eddington-limited accretion onto the black hole. The spherical accretion model (based 
on works by Kazanas, Protheroe and Ellison JO, O) is used in most of the calculations of the neutrino 



production in central regions of AGN |, |l3l yj. According to this scenario, close to the black hole the 
accretion flow becomes spherical and a shock is formed where the ram pressure of the accretion flow is 
balanced by the radiation pressure. We will follow this model even though it is only approximately true 
considering that we will be looking at rotating black holes. 

The distance from the AGN center to the shock, denoted as the shock radius 1Z (~ a few Schwarzschild 



radii) contains the central engine of AGN. The shock radius is parametrized (TT|, [12| by 1Z = xr g where r 
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is the Schwarzschild radius of the black hole, and x is estimated to be in the range of 10 to 100 || (which 
is consistent with the available data ||15|| ). 

The matter density at the shock p(TZ) can be estimated from the accretion rate needed to support black 
hole luminosity, and from the radius and accretion velocity at the shock |J 

p {TV) ~ 1.4 x 10 33 ari LagnQ" 1 gm/cm 3 (1) 

where Q{x) = 1 — O.lx 31 is the efficiency for converting accretion power into accelerated particles at the 
shock [fT2| , and Lagn is the AGN continuum luminosity in units of ergs/sec. In this model the matter 
density falls off as 

_ 3 

Pl (r) = p(K) 2 ; r>1Z. (2) 

Acceleration of protons is assumed to occur by the first order diffusive Fermi jxtechanism at the shock 
E3]. Energy losses during acceleration occur through pp collisions in the gas H; and also through the 



Left-right transitions can occur as a gravitationally induced coherent precession, as was investigated by Cai and Papini, 
, and is in general quite small. 

2 Accelerated protons interact with protons of the accreting plasma to form charged pions. 
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PI ~~ * P + e+ + e ~ an d P7 — > A% processes in the dense radiation fields present, mainly, in the central 
region. All these reactions give rise to high energy neutrinos through the 7r ± — ^ /x^ — decay chain 
[|, [11], HI; H11- These neutrinos are expected to dominate the neutrino sky |§ at energies of 1 TeV and 



beyond. 

Since neutrons are not confined by the magnetic field, they are free to escape from the central core 
M. But neutrons can also produce high energy neutrinos via np and wy collisions. The conclusion of ref. 
M was that neutrons with energy less than approximately 10 3 TeV do escape the central core. On the 



other hand, Stecker et. al. |T9| , found that secondary neutrons will not in general escape the shock radius, 
and argue that a significant amount of power is generated through nj interactions. However, according to 
them, confinement of nucleons in the central core occurs only for energies beyond ~ 10 4 TeV. In our study, 
as we shall see later, the matter density is irrelevant whatever the status of neutrons, confined or not. 

Magnetic effects around AGN have important astrophysical consequences. To estimate a value of the 
magnetic field of AGN which we will use in our calculations, we look at two specific models. Both these 
models assume the "equipartition" condition that the external pressure is matched by the magnetic one. 
According to Begelman et.al. p0|, a characteristic magnetic field B is given by 
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B = m G (3) 

(M/M ) 1/2 

where M is the mass of the black hole considered. Szabo and Protheroe || estimate the magnetic field at 
a shock formed by accretion to be 

B ~ 5.5 x 10 27 Q~5 X -!l A( | n G (4) 

where the symbols have the same meaning, as in equation ([!]). Since the shock is the site of high energy 
neutrino production a good estimate of the value of magnetic field would be that evaluated at the shock. 
In the absence of a detailed model we will assume that the magnetic field remains approximately constant 
inside a region of size equal to the pressure scale height, the average fields of two of these regions will 
be uncorrelated. For L AGN = 10 45 ergs/sec, M = 1O 8 M and x = 10 both of the above expressions give 
B ~ 10 4 G. We will be interested in the vicinity of the horizon where the pressure scale height is ~ r g . 
In this region we will assume that B remains at the above value. Outside the inner core of the AGN the 
magnetic field will drop as (pressure) 1 / 2 . 



3. Dirac Equation in Curved Space-Time and Neutrino Oscillations 

We will study the possibility of neutrino oscillations between different flavors and spins in the environ- 
ment of a strong gravitational field. In this paper we consider, for simplicity, the case of two family flavors 
only. We will also restrict ourselves to the case of Dirac neutrinos interacting minimally (as described by 
the equivalence principle) with the gravitational field. 

The outline of our method to compute neutrino oscillations is the following. Starting with the Dirac 
equation in the presence of a strong gravitational field, we use a semiclassical expansion to different orders 
in the relevant parameters to determine the effective Hamiltonian for positive momentum states. We 
look in particular for the lowest order off-diagonal terms that could possibly induce gravitational neutrino 
oscillations, including the effects of magnetic fields associated with AGN. Once the Hamiltonian is obtained, 
we analyze the possibility of resonances depending on the neutrino mass, energy, angular momentum and 
the black hole mass and angular momentum. The different neutrino transitions and expected neutrino 
fluxes are then investigated. 



3.1 Semiclassical Approximation 

We start with the Dirac equation in curved space-time pifl , 

[<l a {^ + Up) ~ m+ /P L ]?P = (5) 

where are the tetrads, m is the mass matrix, / is the weak interaction current matrix, Pi the left-handed 
projection operator, and the spin connection is 

w„ = IbaM^e" (6) 
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where the semicolon denotes a covariant derivative; we use Greek indices for coordinates in the general 
frame and Roman indices for the local inertial frame. As in flat space, there are four neutrino states for 
each flavor that we separate into two states corresponding to neutrinos traveling towards or away from the 
observer. We call these positive and negative momentum states respectively. 

The first step in the semiclassical approximation is to replace the neutrino spinor by 

i/> = e lS X (7) 
where S is the classical action defined in terms of the Lagrangian density C(x, x) and affine parameter I as 

/ ] d% 
dl£(x,x) C = --g liv x ll x v \ x = — (8) 

(an over-dot will always denote differentiation with respect to I), and solves the Hamilton- Jacobi equation 

g^dpSdvS = 0. (9) 

X is constructed as the solution of the Hamilton- Jacobi equation and determines the trajectory of classical 
particles. 

In order to proceed further we use the approach described by Sakita and Tsani El. If x be the solution 



to the classical geodesic equation for massless particles, that is the classical trajectory, a new set of local 
coordinates 

{Z,£ A }, A = 1,2, 3 (10) 

is chosen such that 

x» = x»(l) + ^ A (l)t A . (11) 

Here solves the variation equation 

*{(£H(©}-° < i2 > 

to 0(£ 2 ). The functions x + vj£ A solve 5S = to second order in £ and so describe (approximately) a 
bundle of neighboring geodesies parametrized by the £ A ; the values of the £ determine the separation of a 
given geodesic in the bundle from the reference geodesic x M (Z). 
The semiclassical momentum is defined as 

Vii = -d^S\ x=s (13) 
where the derivative is with respect to the end-point. The solution to @ is, to 0(£ 2 ), 



s = - P y A i A - \ {K P p»v> p B + \nab)z a z b < 1 1 > 



where T^ p are the Christoffel symbols evaluated at x = x, and 

Nab = ^a u b9^u- (15) 
From the geodesic equation obeyed by it follows that (see Appendix A) 

V^a = C A = constant. (16) 

In the following it will prove convenient to define a time- like vector p which is the component of momentum 
p orthogonal to the namely 

p»=p»-c A {N- l ) AB v» B . (17) 



4 



3.2 Effective Hamiltonian 



With the above preliminaries we are in a position to evaluate the effective Hamiltonian for the neutrino 
system. In order to make a systematic expansion let us define R to be the scale of the metric, so that 
for example, uj^ ~ 1/R] and let p be the order of magnitude of the momentum of the neutrinos. We 
imagine a localized wave packet of extension £ propagating freely through a region of size R, where the 
gravitational field is essentially constant. The uncertainty relation requires that a change in the momentum 
is related to R via Ap ~ 1/R. This in turn implies that the angular spread of the wave packet is given 
by AO = Ap/p ~ l/(pR), so that in a distance R a wave packet spreads RAO ~ 1/p. But this spatial 
spread is determined by the values of £ for which the states are significantly different from zero, so that 



Going back to the Dirac equation in curved space fl5|), after substituting (^), we obtain the equation 

(-e^ a d,S + < 7 a ^ + V -m) X = (18) 

where 

V = *7 a e!H+ !Pl (19) 

and ujp is defined in (|6]). 

We make a double expansion of \i first i n powers of £ and then in powers of 1/pR, as follows 

x = x (o) + xiv + o (a 

x (o) = Uq + Ui + Ui + ... ] U v ~(pR)-» z/ = 0,±,l,--- (20) 



' 2' 

Xa = V 1A + Vs A + V 2A + ---; V uA ~ R-\pRy~» z/ = l,|,2, 



We substitute these expressions into the Dirac equation and make a double Taylor expansion in £ and 
1 / (pR) and demand that each term vanish separately. In the perturbation we choose the mass to be 
(p/R) 1/2 . This is for convenience in defining the perturbation expansion in an effort to look for the 
lowest order off-diagonal terms that might cause gravitationally induced helicity flips of neutrinos. For the 
lowest possible magnitude of the mass, i.e. m ~ 1/R, there are no off diagonal terms in the Hamiltonian. 

The result is the following set of mixed equations for the spinors Uq, Ui • ■ -. To 



p'Uo = 0, p'Ui = mU 



pV x + ^-U + i J A V XA + VqUq -m(Jx=0 
? 

fVz + -Ui + i J A V* A + V Ui - mUt = (21) 

2 -U 2 2 2 

P 

where 

i 

X A = {N- X ) AB {v,b - Cj ^)Kl a - (22) 

p 2 

(the over-bar represents variables evaluated on the geodesic: x — x). To 0(£), 



where 



?^ix=iXA $Vz_ A =UU± (23) 



A a = ? B - e%„P^ A l a . (24) 



From (0) and (|23|) we obtain 



iC/ = 0Ui 



o 



iUr =OUi-— J A V lA - -VoC/o- (25) 

2 2 2 2 
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where 

° = - 1 2^ + ^-^+^J A U (26) 

It is now possible to reduce these equations to a Schrodinger-like equation (as shown in Appendix B) 
involving only \ which reads 

O - - V J x- (27) 

The above equation describes a set of states which are almost pairwise degenerate (due to the condition 

f>X — 0); the effective Hamiltonian H e g for the positive momentum states is obtained using degenerate 
perturbation theory. Denoting by P + the projector onto the positive momentum states we find 

H eS P + = P + (O--V )P+. (28) 

It is always possible to go to the p rest frame where p a = E(l, 0,0,0) and p a = E(l, 0,0,1), neglecting 

the neutrino rest mass and using p 2 = p ■ p; in this (local Lorentz) frame, the projector onto positive 
momentum states is 

P+ = diag(l, 0, 0, 1) (29) 

In what follows we choose frames where the momentum has the form given above. In Appendix C we show 
that, with these previous considerations, equations (^7]) and (|28|) lead to 

Hcs = ia + ^m 2 -p- J cS P L + mQ b r b (30) 

where a is spin and flavor diagonal and has no observable consequences i. The effective current is 

JcS = J\v ~ Jq- (31) 
Here is the weak interaction current, Jq is defined as 

J a c = \^% cd ^ + (32) 

and 

A/ C d = (e /AM , - e/^e^, r) ah = diag(l, -1, -1, -1) (33) 

In order to determine we note that in the rest frame of the accreting matter it takes the form 
Jw;v e = Gf{^N p — N n )/y/2 and J^-,^ = —GpN n /\/2 where Gp is the Fermi coupling constant and 
N p and N n are the proton and neutron number densities respectively (only the zero-th component is 
non- vanishing in this frame). In the neutrino frame used above this becomes 

r w . ve = ^2N P -N n )u a , j^ = -^£ Nn u a (34) 

where u a denotes the four velocity of matter in the chosen local inertial frame. In the model considered, 
matter is ultrarelativistic (within the region of interest), so that u 2 ~ 0. 
Finally, 

e' = \V P -f + \ (JS. + \^) U + w> '^"' ) ■ P5) 

In addition to the gravitational effects it is essential to incorporate the effects of the magnetic field 
associated with AGN. It is well known that if a neutrino has non- vanishing magnetic moment (or transition 



3 An imaginary contribution to a would signal a decrease in the overall neutrino flux and is also unobservable. 
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magnetic moment of dipole type) its helicity can be flipped when it passes through a region with a 
magnetic field which has a component perpendicular to the direction of motion. This idea has been 
analyzed by Okun et.al. ||24|| , and more recently the combined effects of flavor mixing, magnetic spin flip 

H, El- 



and matter interactions have been considered 
stems from a term of the form 



The interaction with the electromagnetic field 

fia ab F ab iP (36) 

(F ab is the electromagnetic field tensor and a ab = |[7 a ,7 6 ]-) to be added to the left-hand side of (|5|); just 
as in the above references only the electric and magnetic fields in the direction orthogonal to p contribute. 
Denoting by E r and B r the components of the electric and magnetic fields respectively in the frame where 
p a = E(l, 0, 0, 1), and using (|3~6|) we find that the effective Hamiltonian is modified according to 



H. 



cff 



H eS + HVP-P 




n 



(E 1 + B' 



(37) 



where /i denotes the magnetic moment matrix and B a , E a denote the magnetic and electric fields measured 
by a locally inertial observer. 

We consider two generations with four-component Dirac neutinos. Here we examine the v e — system. 
The same can be done for the v e — v T and — v r systems. Using the chiral bases u eL , v^ L , u eR , v flR1 we 
can write the evolution equation for neutrino propagation through matter in the presence of a strong 
gravitational field as H 



.d 



I 




\ 




f 




\ 




























V 




J 




V 







(38) 



where H e Q is the 4x4 matrix (0), containing the effects of the weak and gravitational interactions and 
the effects of electromagnetic field. An order of magnitude estimate (supported by explicit computations 
for particular metrics) show that the terms (eq. |35], off-diagonal terms due to gravitational effects) in 
the effective Hamiltonian are negligible compared to the effects of the magnetic field for all interesting 
cases (see section 3.4). The gravitational effects in the diagonal elements however are very relevant. The 
final expression for the effective Hamiltonian H c q is the following 



/ -p ■ Jc| + |Am? 2 sin 2 $ 



H, 



cff 



|Amj 2 sin 2$ 



Efj, ee Q, 

En^Vt 



|Am 2 2 sin2^ Efi ee Q* Efj^Q* \ 



-p- J e ^ + |Am 2 2 cos 2 ^ Efi^n* Efi^Vl* 

E/i etl n o o 

%^ iAm 2 2 J 



(39) 



where d is the neutrino mixing angle, J e g is given in (|3l|), E is the energy of the particle, and Q is defined 
in (p7|). Note that, because v eR and u m are sterile, i.e., do not interact electroweakly with matter, they 
can be considered vacuum mass eigenstates. 

3.3 Resonances 

Using (|39|) we can determine the AGN regions where resonant transitions occur (a phenomenon remi- 
niscent of the MS W effect @, ) . These resonances are governed by the 2x2 submatrices of (|39D for each 

4 Okun p4j] pointed out that the manifestations of an electric dipole moment and magnetic dipole moment are practically 
indistinguishable for a highly relativistic neutrino. For simplicity the term "magnetic moment" will be used in this paper to 
represent the combined effect of an electric and magnetic dipole moment. 

5 Note that the left hand side involves differentiation with respect to the affine parameter which has units of E~ 2 the units 
of length being E^ 1 . Therefore H e s has units of E 2 which differs from the usual units of the Hamiltonian. 
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pair of states B. 

The five possible resonances are obtained by equating the diagonal terms for each submatrix and give 
rise to the following resonance conditions 

^Am 2 2 cos2tf + \pl (p ■ u) GpNp = (z^-^J (40) 

^Am 2 12 cosH + &-^-G F N n + p-J G = (zv-^O (41) 

-iAm? 2 sin 2 ^ + VL^-G F N n + p ■ J G = iy^v^) (42) 

iAm^cos 2 ^ + G F {2N P - N n ) - p ■ J G = {y eh -+v m ) (43) 

-iAm 2 2 sin 2 tf + ^5 G F (2N P - N n ) - p ■ J G = {y e ^v eR ) (44) 

As will be shown in the following section the matter effects are very small compared to the gravitational 
effects which makes J G dominant whenever present (this makes irrelevant the precise value of N n for the 
cases (41-44); only for the MSW resonance 0, ||, eq. fl40f) is the matter density important). 

In order to determine whether the above resonances induce significant transition probabilities we con- 
sider the corresponding 2x2 submatrices in detail. Each of these can, by a suitable subtraction from the 
diagonal term, be cast into the form 

(t! d ) <«> 

where d is related to the matter and gravitational part and b to the magnetic field; their explicit form is 

b = Efi„n d = -p ■ J cS + (Am 2 term) (46) 



where E, \x v and Q are respectively the energy, magnetic moment of the neutrino and electromagnetic field, 
see (PTD; the "Am 2 term" depends on the specific 2x2 matrix and can be easily obtained from (p5|); it is 
always smaller than |Amf 2 |- 

Resonances occur when d vanishes, in which case the transition probability is well described (for suffi- 
ciently slowly varying d and b) by the Landau-Zener approximation ^Z7\ 



^ = exp^-27r 2 ^) (47) 



where 

(3 = b | r es a = d | res (48) 

The condition for these resonances to induce an appreciable transition probability is 

I? > £. (49) 

The presence of a magnetic field can also induce coherent precession of the states which has been 
studied in sufficient detail in The condition for this to occur is, for slowly varying b and d in (fJSI), 



V ° £ R > 2tt (50) 

where R is the magnitude of the region where the magnetic field is coherent. In this case the transition 
probability is 

1 b 2 

Pprcc " 2~WTd? (51) 
which is significant when b > d. So, while the presence of a diagonal term requires a shorter distance 
R to generate rapidly varying phases, the very same effect decreases the transition probability: neutrino 
interactions with matter or gravitational effects effectively quench spin precession |23[ . 

6 We will ignore multiple resonances which occur whenever three (or four) diagonal elements concide; such resonances can 
occur only under very restricted circumstances. 
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3.4 Estimates 



The orders of magnitude for the weak and gravitational currents are, from ([51) and 

J w = "fZJL „ i - 33 p eV' 1 , J G ~ R- 1 (52) 
m p 

where p is in units of eV 4 . According to ([!]) and (J^) the order of p for typical cases is 10 1 — 10 4 eV 4 . Clearly 
from above and taking R ~ r g , the gravitational current part is found to dominate the weak current part 
for all relevant values of r g (10 14 to 10 20 eV _1 ). 

The order of magnitude for the m6 term in the Hamiltonian will be (considering that ~ 

mQ b ~ mR~ x ~ mr^ 1 . (53) 

If we compare this term with Ep u Q for typical values (E ~ ITeV, Q ~ 10 4 G, r g ~ 10 18 eV -1 ), for a very 
small magnetic moment p ~ lO -19 ^ and neutrino mass m ~ 1 eV, we get 

E ^ 10 6 . (54) 



Therefore, in what follows we will neglect the term. 

The values for |Am 2 2 | corresponding to different scenarios of neutrino oscillation are 





Am 2 2 cos 2 $ (eV 2 ) 


Am 2 2 sin^ (eV 2 ) 


Vacuum 




l(T iu 


io- n 


Solar small angle 


29] 


lir 6 


10- 8 


Solar large angle 


J29|| 


lir 6 




At. neutrinos |29 




10- 3 


10- 4 


lsnd g 


1 


lO" 3 



Table 1: Approximate values for Am 2 terms. 



The MSW resonances ( KDj ) correspond to 

Am 2 2 ~ l(T 33 £p (55) 

which, since p < 10 4 eV 4 and E < 10 19 eV (in order to get an appreciable flux) corresponds to Am 2 < 

lO -10 eV 2 . Thus, aside from these extremely small values, the usual MSW scenario does not take place. 
Since p ~ E, neglecting weak current part, 

E 

V JeE ~ — (56) 

For a typical value of energy E ~ 1 TeV, beyond which AGN neutrinos start dominating the sky, and 
r g ~ 10 20 - 10 14 eV" 1 (M = 1O 1O M - 1O 4 M ), p ■ J eff ~ 10" 8 - l(T 2 eV 2 . The orders of magnitude 
from the previous table show that values of Am 2 corresponding to the solar and atmospheric entries would 
undergo resonances in the vicinity of the black hole. Larger values of the mass difference (such as the 
LSND entry in the above table) would resonate only at significantly larger energies (above 100 TeV); in 
order to have comparable resonance for Atmospheric and LSND |Kj values we need even more energetic 
neutrinos (~ 10 2 TeV and ~ 10 6 TeV respectively). A point to note here is that all resonant transitions 
do not occur simultaneously as can be seen from equations (H0-H4T). At fixed conditions for instance, the 



transitions (v Mi — >v eR ) and (v eL —*u eR ) cannot both be realized together. 

In order to determine whether the above resonances induce large transition amplitudes we will need to 
estimate the magnitude of Q. The models we consider all assume that the magnetic field is determined 
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by the "equipartition" condition M that the external pressure is matched by the magnetic one. The scale 
height for the pressure is ~ r g and so we expect the magnetic field to be uniform only through distances 
of order r g ; on larger scales it will vary randomly. As we mentioned in section 1, the value of the magnetic 
field within these "coherence patches" is \Q\ ~ \B\ ~ 10 4 G. Let us now define the scale height A for the 
effective current terms 

dl 

A 



d(\np ■ J cfr ) 

which has units of the affine parameter I. Then a ~ Am 2 /A and the condition ( T4D| ) becomes 

1/2 



(57) 



1 



EB 



Am 



2tt 2 A 



Resonant transitions will occur then if the relevant magnetic moment (or transition magnetic moment) 
satisfies this bound; this constraint depends, through A, on the type of metric and will be studied for some 
cases of interest in section 4 below. Note that we have assumed that the magnetic field remains constant 
over an interval of magnitude ~ A. 

Within this estimate the condition for coherent precessions fl5"0|) becomes 



(pBr g ) 2 + (*?pL + If > 4vr 2 (59) 



and the transition probability Q5"T| ) becomes large provided 

1 / Am 2 r 



"" - B7 g (J + ^PJ = ^ (60) 



which is the condition for coherent precession to generate a significant number of helicity flips. 
Comparing the different magnetic moment expressions obtained above we get 

V2 / A •> \ -1/2 



( Am 2 r g \ ( Am 2 r 



<:l 



/^min V E J V E 



(61) 



which implies /i^~2/i^; n for any value of Am 2 r g /E. 

Though realistic calculations would require an accurate profile for the magnetic field, the above results 
do provide useful order-of-magnitude estimates of the relevant quantities. We will use these expressions in 
section 4.2 to estimate the possibility of transitions of each type. 

4 Neutrino oscillations in a Kerr space time. 

In this section we will apply the formalism developed to the interesting case of a Kerr black hole. 

4.1. Hamiltonian for the Kerr metric. 

The gravitational field of the rotating black hole is given by the following axially symmetric stationary 



Kerr metric f31 |: 



1 2rgra sin 2 6d(f) dt (62) 



P 2 



where 



A = r — r g r + a , p = r + a cos 8 (63) 
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This metric depends on two constant parameters M and a, where M is the mass of the black hole and 
a is related to the angular momentum by the relation Lbh = ma; r g is given by r g = 2M in geometrized 
units (c = G = 1). For a = the Kerr metric becomes the Schwarzschild metric in the standard form. 

We now consider the motion of a particle of mass m in a Kerr field 0. The Hamilton- Jacobi equation 
is fully separable in this case [^]. Writing the action S as 

S=-Et + L4> + S r (r) + Sg(9), (64) 

where E is the energy and L the angular momentum of the particle. We obtain two ordinary differential 
equations 

2 i 2 

+ (aE sin 9 -) + a 2 m 2 cos 2 9 = K (65) 

sin 9 

- ^{r 2 E + a 2 E - ahf + m 2 r 2 = -K (66) 

where K is a new constant of motion. The four-momentum of the particle is then 

a 1-I77 -ic/ -lc/ L — g 03 E/g 00 \ 

P = I 7o E , 7i 5 r , 72 5 e , I (67) 

where 

7o = v^oo 7i = V~9ii k 2 = — ~ 933 (68) 

#33 

Our calculations involved a set of tetrads in which the momentum takes the form p a = E(l, 0, 0, 1). The 
corresponding tetrads are 

e a 

A' 

where 

s a = sin a, c a = cosa (69) 

and similarly for /3; a and f3 are the polar and azimuthal angles of p a in (|57D. The angles fulfill the following 
conditions 

1 2 

tana = — - tan/3 = . (70) 

P 3 ^( p l)2 + (j9 3)2 

Using the above relations we get an expression for the effective Hamiltonian for Kerr metric (|30|). We 
use this to look at resonances in the next section. The explicit expression opf the various terms in the 
Hamiltonian is quite involved, we include a brief description of such terms in Appendix D; here we present 
the result of numerical analysis of the expressions. 

As a limiting case of the previous analysis, we consider the limit a — > 0, corresponding to a Schwarzschild 
black-hole. The neutrino geodesies lie on a plane which we take as the 9 = 7r/2 plane, we then obtain 
P u Jq = and the effective Hamiltonian reduces to 

1 2 EJw° /0 0\ / Q* \ 
H >—2 m2 -J=77r[0 l) + "U ) <71) 

where we assumed the absence of matter currents and where m and /i denote the mass and magnetic 



/ 7o 








V 












JvSq, 









72C/3 


kspc a 




V o 




-72.S/3 


kcpc a 


) 



moment matrices respectively. Thus, except for a factor 1/yl — r g /r which is important only for r ~ r g 
the above expressions reduce to the flat-space situation. In particular there is no gravitational contribution 
to the diagonal elements. One can also easily determine the off-diagonal terms Q b in (0), which give a 
contribution mLr 2 /(2.Er 2 ) to H e ^ (r 2 denotes the usual Pauli matrix). As mentioned previously these 
terms are negligible for all cases of interest due to the factor of m/E. 



7 Included only for completeness and to define the notation used later. 
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4.2 Resonances for the Kerr metric 



We need to determine whether resonances are possible in the vicinity of a Kerr black hole. For a rotating 
black hole the spin flavor transition is due to the transference of orbital angular momentum to spin angular 
momentum and also due to the transference of the black hole's angular momentum to the neutrino spin. 
To determine possible resonant behavior, we will find the regions (r, 9) where each of equations (|H) to 



J) are satisfied as a function of the parameters E, L, K, a, r g . We will do this for different values of AmJ 2 
and mixing angle. We choose normalized parameters 



J 



Er n 



K 



{Erg) 2 



(72) 



for a neutrino obeying the equations (|65[) and fl66|). It is easy to see that p ■ J eS can be written (neglecting 
the matter terms) in terms of the energy and Schwarzschild radius in the form 



\p- Jeff| = Er g l f(r/r g ,6,j,k,a/r g ). 



(73) 



In Fig. 1 we show In |/| as a function of (r, 6), for an allowed value of the pair (j, k), where an angular 
momentum (a = 0.4r 9 ) has been chosen for the black hole.Q We have plotted the magnitude of p ■ J e g in 
an effort to give a clear picture of how resonances take place in this particular metric. It is to be borne in 
mind however that all resonance transitions do not occur simultaneously. Figure 2 presents the / contour 
plots under the same conditions as in Fig. 1. 




Fig.l. Plot of In |/|, see @, for j = 0.4, k 

1O 8 M . 



0.15, a/r g = 0.4 and M 



8 A realistic analysis of accretion onto black holes must account for the fact that the central hole is quite probably a Kerr 
black hole with angular momentum parameter a only slightly less than the gravitational mass p3] . 
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e 



Fig.2. Contour-plot of /, see @, for j = 0.4, k = 0.15, a/r g = 0.4 for the 
case M = 1O 8 -M0. The contours a,b,c,d,e correspond respectively to /(g) = 
±1, f{b) = ilO" 1 , f(c) = ±10~ 2 , /(d) = ±10" 3 , /(e) = ±10- 4 , where 
positive values take place for < < 7T/2, and the negatives ones for 7r/2 < 

< 7T. 




Figure 3 illustrates the dependence of the function / for a specific value of (r, 9), on allowed values of 
(j, k) and taking a = 0Ar g . 

Resonances occur provided / is comparable to ±Am^ 2 cos 2 $r g /E or ±Am 2 2 sin 2 $r g /E as can be seen 
from equations ( (41~|) to (j44|) . As an example we can analyze the case of energy E ~ 1 TeV, r 9 ~ 10 18 eV _1 
and consider the Solar large angle solution, Am 2 2 ~ 10~ 6 eV 2 . Comparing these with / from figure 2, we 
find that it corresponds to contour "a" . We conclude then that resonances occur in the vicinity of the AGN 
for this choice of energies. In the same way it can be seen that resonances are present for all the given 
values of Am 2 in Table 1, all relevant black- hole masses and angular momenta, and all neutrino energies 
above ITeV. 

It is also of interest to determine the energy dependence of the resonance conditions, which is presented 
in figures 4 and 5. 
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20 



Fig.4. Plot of log |/|, see (|73|), as a function of T and E for the case 9 = 7r/4, 

j = 0.4, fc = 0.15, a/r g = 0.4 and M = 1O 8 M . 



loglf I 




E (TeV) 



Fig. 5. Plot of log | / 1 , see @, as a function of 6 and E for the case T 

j = 0.4, k = 0.15 and a/r g = 0.4 and M = 10 8 M G 



6r n 



For fixed values of j, k, r and 6, the function / depends linearly on the energy E. 
In order to determine the values of the magnetic moment which will allow for these resonances to induce 



large transition probabilities, we defined the quantity /i^ n in equation 



This is plotted in Figs. 6 and 



7 for some representative values of Am and E and for a magnetic field of 10 G. 
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Fig. 6. Plot of 10 13 /lmin> see (§|), for the case Am 2 = eV 2 , E = lTeV, 

j = 0.4, k = 0.15, a/r g = 0.4 and M = 1O 8 M . 

The resonances which we described above will induce large transitions provided //„ is larger than //^f n ; 
for the example considered this corresponds to ~ 10~ 13 — 10 _14 /ie (where \ib denotes the Bohr magneton). 
The above requirements lie comfortably inside the direct experimental bounds (fi u < 10~ 10 jjb) [HI as well 
as the astrophysical limits {fi v < 10~ u [1b) [0 and such values for the magnetic moment are consistent 
with a wide variety of models ||35|| . In view of this the above resonances will induce significant transition 
probabilities whenever the resonance conditions ( f40| - fi4l ) are satisfied. 




Fig. 7. Plot of 10 13 /i™- n , see (H), for the case Am 2 = 10~ 6 eV 2 , E = lTeV, T = 

6r g , 9 = 7r/2, a/r g = 0.4 and M = 1O 8 M . 

The minimum of / /i^L (HI) (note that equation ( |6l]) involves only an order of magnitude calculation 
where |/| ~ 1, (|73|)) corresponds to Am 2 r g /E ~ 1, and then /i^fn is comparable to //^f n ; for E = ITeV 
and r g ~ 10 18 eV _1 this corresponds to Am 2 ~ 10 _6 eV 2 which are the parameters used in plots 6 and 
7. Explicit calculations show that > /^in whenever |/| > 1, which at resonance corresponds to 
Am 2 r g /E > 1. The transition efficiency is greatest for the minimum value of the magnetic moment and 
corresponds to the spin flip in vacuum. 

In order to estimate the effects of these resonances on the neutrino fluxes on Earth we note that, 
neutrinos are expected to be generated mainly through the n — » ji — > e decay chain, implying that we 
can expect twice as many muon-type neutrinos as electron-type neutrinos. A negligible number of tau 
neutrinos are produced in an AGN environment. Even though the above calculations were done for the 
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case of two species, we expect the same results for the case of three species: resonances will occur for all 
possible transitions. Such resonances will tend to "spill over" the neutrinos from the highly populated 
states into the less populated ones evening out the neutrino numbers in all states. It follows that, with 
these simplistic assumptions, the muon-neutrino flux will be decreased by a factor of 1/4, the electron 
neutrino flux will decrease by a factor of 1/2 but we can expect as many tau-neutrinos. 



5. Conclusions 

The evolution of neutrinos in an AGN environment is susceptible to gravity-induced resonances. Such 
resonances would induce various spin and spin-flavor transitions whenever equations Q4"TI-fP|) are satisfied. 
We have shown that such conditions would be expected to be satisifed for the currently accepted AGN 
models and for various neutrino parameters currently used in the literature. 

The lack of precise modeling of the AGN magnetic field (eg. is the field dispersed or is it concentrated 
in flux tubes), and lacking a better understanding of the parameters in the neutrino system it is impossible 
to determine unambiguously whether such resonances would in fact occur in a given AGN. 

Using the estimates for the magnetic fields, and approximating the gravitational field by the Kerr 
metric, the various values for Am 2 we use do lead to resonances of the type ([Tj]- |44]) . The corresponding 
transition probabilities will be large for a sufficiently large magnetic moment. Each resonance will occur at 
a different place along the neutrino trajectory; plots 1, 2 and 3 show that essentially all neutrinos produced 
near the AGN core will traverse regions where all four resonance conditions are satisfied (for the neutrino 
parameters under consideration). Such resonances induce large transition probabilities for a wide variety 
of neutrino masses and energies (provided the neutrino magnetic moment is sufficiently large). 

We do not expect resonant transitions to be 100% effective, whence the neutrinos will tend to populate 
all available states. Lacking detailed knowledge of the AGN environment, but relying on the fact that all 
resonances are viable, we can only estimate that all such states would end up by being equally populated. 
This, of course, depends on the masses of the neutrino states. For very small values of Am 2 , Am 2 < 
10- 10 eV 2 , © will be satisfed and MSW resonances would occur; the other posibilities (|41|-444]) can be 
realized, but now farther from the AGN core where the estimate of the gravitational field is less reliable 
and where matter effects can no longer be ignored. We will not consider this (very complex) possibility 
further in this paper 

It is interesting to note that for Am 2 > 10 _10 eV 2 , the resonances described in this paper are not 
related to the MSW resonant flavor oscillations (which are suppressed in this environment due to the low 
densities present unless Am 2 < 10~ 10 eV 2 , see sect. 3.2), but to the interplay of the gravitational and 
electromagnetic interactions with the mass terms. Also it must be noted that all resonance transitions do 



not occur simultaneously as can be seen studying equations (f|0|-[44|). 

Noting that the AGN neutrinos travel through cosmological distances before being observed on Earth 
and taking the matter effects to be negligible on the way to the earth one might argue that the vacuum 
oscillations generated by the non-zero mass matrix would produce identical effects to the ones considered 
here. For a distance of approximately 100 Mpc (a sensible distance to AGN), vacuum oscillations occur for 
Am 2 < 10" 19 eV 2 for 1 TeV energy neutrinos which is restrictively small. The effects of a cosmic magnetic 



field (the typical value of galactic [3~7| and galactic cluster |38|] magnetic fields is of the order of 10 6 G; the 



intercluster magnetic field is <10~ 9 G [|39| |) can, of course, induce rapidly varying phases in the neutrino 
states, but the transition probability (|51|) is nonetheless small. 

In our calculations we have chosen a rapidly rotating black hole (a = 0Ar g ). This assumption is based 
on the conclusions of reference |[33|| . For lower values of a the conditions under which the resonances 
discussed in ths paper can occur become more stringent. For the extreme limit of a Schwarzschild back 
hole, transitions are expected only in the immediate vicinity of the horizon. 

In conclusion we can say that gravitational induced resonance transitions are the most important pro- 
cesses which can generate effective spin flavor conversion in Active Galactic Nuclei. The precise expressions 
for the flux to be detected at Earth-bound neutrino telescopes depends, of course, on the details of the 
AGN model. We have shown, nonetheless that for a very wide range of neutrino masses and even with 
very small magnetic moments, neutrinos will undergo resonant transitions. In their trip through the AGN 
environment the neutrinos can experience several such transitions which tends to even-out the population 
of all neutrino states. On one hand this decreases the expected flux of electron and muon neutrinos, on 
the other hand it dramatically increases the number of tau neutrinos of energies of 1 TeV and above. 
This provides strong motivation to search for such tau neutrinos by the existing neutrino telescopes, and 
the recently proposed 1 KM3 detector |j3~6|| . Even if the magnetic moments are ~ 10 _14 /^b we can expect 
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an observable number of, for example, 10 TeV r neutrino events; the precise number of such events is 
determined by the flux of electron neutrinos at that energy. 
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Appendix A 



In this appendix we give a derivation of equation fllED, Consider a geodesic x^{l) for which the classical 
momentum is x^ = p^ and solves the Hamilton- Jacobi equation g^p^p 1 " = 0. Taking an I derivative we 
then obtain 



. dg a/ 3 g 
2P » = ~dx^ P P 



(A.l) 



Expanding now x to first order around the geodesic, x^ = x^ + (recall that this expression describes 
geodesies to 0(£ 2 )), and substituting in the Hamilton- Jacobi equation yields , 



2p,e + ^p a p"^ = °' 

where £ M = v\^ A . From (|A.1|) and ( |A.2|) it immediately follows that 

d 



dl 



(A.2) 



(A.3) 



Using the fact that the £ A are independent of I it follows that p^v^ = ca =constant as was to be shown. 



Appendix B 

In this appendix we show how to obtain a Schrodinger-like equation from equation Q2"5|). We note that 



But, from (GO 



iU + iUi = O(U + Ui)-— /C Via - -V U . (B.l) 



iU + iUi = ix + 0(l/R 2 ) 

yT% = ^JC A d AX + 0(l/R 2 ) (B.2) 



and therefore by reordering terms in( |B.l|) 

i(d t + f X^) X = Ox + 0(1/R 2 ) (B.3) 

where O is given by 

~ 777, _ 

O = O-jV , (B.4) 

The second term in ( |B.3|) describes the spreading of the wave packets and is not relevant for our discussion. 
We therefore ignore such terms. The final equation therefore reads, 



i X = Ox (B.5) 



Appendix C 

In this appendix we evaluate the effective Hamiltonian, given by equation (p8|) , H e gP + = P + (O — mV /2j P + 
where O, Vo are defined in (EH) and (|19l). The explicit expressions are 

77? 2 — 1 

O = ia + — -ip a % bc a bc - Pa (J a --e M % cd )P L 

Vo = (-^+^ a )7a-(^ C S C d+^J a )7a75 (CI) 



20 



where 



= in ,7 ]■ 



(C.2) 



Using the explicit form of P + and working in the chiral representation for the gamma matrices we obtain 



p +1 a p + = ie abcd n ® r 1 

p 2 



P+a ab P+ 



:- v « b + papb+p y b - papb )n®sp, 

p 2 

1 



2p 



(p a p b - pY -i€ abcd p c p d j 5 )P + 



when p a = E(l, 0, 0, 1) and p a = E(l, 0, 0, 0). Using these relations we obtain 

H eS = ia + ^m 2 - p ■ J cS P L + mQ b r b . 
where J e s is given by (|32 ) and Q b by (|35|) . 



(C.3) 



(C.4) 



Appendix D 

In this appendix we show how to calculate the p ■ J term in the Hamiltonian; the metric is given in 
section 4.1. We will use the following notation 



2 _ #03 



903 



I I ; 2 ^U3 3uo 

7o = \/9oo 7i = V -0« fc = #33 V= — 



along with the tetrads 



#33 

/7o 77 \ 

7iC Q A;s Q 

V -7lC/jS a -72S/3 /cC/jCq / 
6 



7o 



and 



/ 7o _1 

-kECaS/3 7r 1 C/3 -72~ 1 S a S / 3 -k~ l C a Sp 



where 



- — 903(900933 ~ 9m) 1 
With these expressions we construct A a ft c according to 

Explicitly 

p. J = p. (J w + J G ), 



;d.i) 



(D.2) 



(D.3) 

(D.4) 
(D.5) 



where the currents are given in ( |32| ) and (|34"D, being the weak current part negligible compared to the 
gravitational part (as shown in subsection 3.4), and 



V ' Jg — 2^(^132 + -^213 + 3A 32 i + 3A 21 + A102 + A 2 10/ 



(D.6) 
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where, for example, 



A 132 = [^(k - c a c pl \)d x + ^-(ksl + 72 2 4)9 2 + ^(1 + kc,)d 3 ]ks p . (D.7) 



The angles fulfill 



1 2 

tana = -^- tan/3 = P (D.8) 



P 3 ' ^/(pT + Cp 3 ) 2 ' 



In the expression for A a fe c , 9j with i = 1, 2, 3 represents d r , d$ and respectively. Similar expressions hold 
for the other values of X a b c , which we will not list. 



22 



